Magnetic phase separation and superfluid density measurements in alkali 
metal-organic solvent intercalated iron selenide superconductor Li x (C 5 H 5 N) ;/ Fe2_ 2 Se2 
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We report the low-temperature electronic and magnetic properties of the alkali metal-organic sol- 
vent intercalated iron selenide superconductor Li a: (C5H5N)j / Fe2-zSe2 using muon-spin-spectroscopy 
measurements. The zero-field /iSR results indicate that nearly half of the sample is magnetically 
ordered and there is a spatial phase separation between the superconducting and the magnetic frac- 
tion. The transverse-field /iSR results show that the temperature dependence of the penetration 
depth A (T) in the mixed state of Li a; (C5H5N)j / Fe2_ z Se2 can be explained using a two-gap s + s- 
wave model with gap values of 6.82(92) and 0.93(7) meV. This implies that the symmetry of the 
superconducting gap in this system remains unaltered to the parent compound FeSe even after the 
intercalation with the molecular spacer layer. We obtain A(0) = 485(21) nm at T = K. 

PACS numbers: 76.75.-K, 74.70.Xa, 74.25.Ha 



Since the discovery of the iron-based superconductor 
LaFeAsOi-zFj, in 2008 [1], it has been fascinating to 
observe how the research evolved during the last four 
years. In general there are two classes of iron-based su- 
perconductors, chalcogenides and pnictides. It was the 
iron pnictide system which dominated for the first couple 
of years due to their higher T c values, with a highest T c 
of 56 K for Gdi-zThcFeAsO 0. Now, the iron chalco- 
genide system has generated lot of interest to the con- 
densed matter physicists with the development of new su- 
perconductors with more and more higher T c values and 
a fascinating coexistence and competition with strongly 
magnetic phases . The parent compound of this sys- 
tem is FeSe, superconducting with a T c of 8.0 K at ambi- 
ent pressure [l(| and 37 K at 7 GPa [Hj]. The T c value can 
be raised as high as 65 K for a single layered FeSe films 



grown on SrTiC>3 substrate [12h14| . The substitution of 
tellurium on the selenium site also increases T c to maxi- 
mum of 14.5 K at ambient pressure and 23.3 K at 
3 GPa [l3| ■ A further milestone in the evolution of iron 
chalcogenide superconductors has been achieved by in- 
tercalating alkali metals (K, Cs, Rb) between FeSe layers 
which increases the T c values above 30 K 

SHE!- These 

alkali intercalated materials also have a second supercon- 
ducting phase beyond 12 GPa with a T c of 48 K 0. 
Interest has been redoubled in this class of superconduc- 
tors with the recent discovery of the enhancement of T c of 
FeSe above 40 K by the intercalation of molecular spacer 
layer 2(3- 23[. It has been claimed that the value of T c 
increases with increasing inter-layer distance between the 
FeSe layers 24|. This has stimulated a tremendous re- 



search interest to elucidate the superconducting pairing 
mechanism, and gap symmetry and how they evolve with 
the increasing inter-layer distances in these intercalated 
FeSe superconductors. 

In this letter, we report detailed investigation of the su- 



perfluid density and the symmetry of the superconduct- 
ing gap in the alkali metal-organic solvent intercalated 
iron selenide superconductor Li a; (C5rl5N), / Fe2-zSe2 by 
muon spin rotation/relaxation (/iSR) technique. Our 
results clearly demonstrate that nearly half of the vol- 
ume fraction of Li. E (C5H5N), / Fe2_zSe2 is in a magnet- 
ically ordered state. The other half is paramagnetic 
and becomes fully superconducting below T c . The ob- 
served A -2 (T) is found to be well described by a similar 
two-gap s + s-wave model, as seen in the parent com- 
pound FeSeo.85 [251 ]. Comparing these results with an- 
other newly discovered intercalated iron selenide super- 
conductor Li .6(NH2)o.2(NH 3 )o.8Fe2Se2 [21|, we suggest 
that the 2 dimentional (D) FeSe layer generate supercon- 
ductivity in this class of materials and the value of T c 
depends only on the superfluid density within this FeSe 
layer which does not change with the increased interlayer 
distance. 

/iSR measurements were performed using the general 
purpose surface (GPS) muon instrument located on the 
7rM3.2 beamline of the Swiss Muon Source at the Paul 
Scherrer Institute, Villigen, Switzerland. For details 
about this technique, see Refs. [2^-31|. Data were col- 



lected both with zero field (ZF) and transverse field (TF) 
modes. In the TF mode, an external magnetic field 
(H = 30 Oe) was applied perpendicular to the initial 
direction of the muon spin polarization. The magnetic 
field was applied above the superconducting transition 
temperature and the sample then cooled to base tem- 
perature, i.e. so-called field-cooled (FC) procedure. A 
continuous-flow 4 He cryostat was used to collect the data 
between 1.6 to 250 K. The number of positron events were 
18 million for each data point. 

A polycrystallinc sample of Li a; (C5H5N)j / Fe2-zSe2 
(x pa 1.0, y w 0.2 and z ps 2.0) was prepared via room 
temperature intercalation into the iron selenide matrix in 
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FIG. 1: (Color online) Temperature dependence of 
the ac-magnetic susceptibility (real part only) of 
Li;c(C5H5N) H Fe2_zSe2 measured using zero-field-cooled 
mode under a magnetic field of H ac = 1 Oe. The inset 
shows the diamagnetic transition in LL(C5H5N) y Fe2- z Se2 
at 40(1) K. 



pyridine solution of the Li metal as described in Ref. [23| . 
A circular disk shaped pellet of the powder sample of 
Liz(C5H5N)j,Fe2- z Se2 with 8 mm diameter and 2 mm 
thickness was mounted on a sample holder. The pellet 
was covered with a thin layer of polymer and kept in an 
inert atmosphere to inhibit air decomposition. 

Temperature dependence of the ac-magnetic suscepti- 
bility of Li x (C5H5N)j,Fe2- z Se2 was performed in a Quan- 
tum Design PPMS magnetometer under a magnetic field 
of H ac = 1 Oe, using zero-field-cooled (ZFC) protocol. 
Figure Q] shows the real part of the ZFC ac susceptibil- 
ity of Li x (C5H5N) y Fe2- z Se2, measured at temperature 
ranging from 5 K to 300 K. The inset of Fig. [T] shows 
a diamagnetic transition in the susceptibility data with 
a T c onset at 40(1) K. A rather large positive signal in 
the normal state susceptibility data shows the presence 
of a sizeable amount of ferromagnetic impurity in the 
sample. This may be due to the presence of some small 
clusters of Fe ions and other Fe-based magnetic phases 
such as a-FeSe, FerSeg, Fe3C>4 in the sample. Comparing 
the normal state susceptibility value with the saturation 
magnetization of elemental Fe, Fe7Ses and Fe3C>4 we es- 
timate that the amount of this magnetic impurity is not 
more than 1 % of the sample mass. A similar scenario 
has also been observed in many other Fe-based supercon- 
ductors 
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In order to search for any magnetic anomalies (static 
or fluctuating) in Li x (C5H5N)j,Fe2- z Se2, we have per- 
formed a ZF-/iSR study. Figure [2] (a) shows the ZF- 
/iSR signals collected at 1.6 K, 50 K and 250 K for 
Li a; (C5H5N)j ; Fe2- z Se2. The signals are almost identical 
in nature and there is no precessional signal and no big 
change in the observed relaxation rate between data col- 
lected above and below T c . However, we have observed a 
large drop of the asymmetry value of the ZF /iSR signal, 
occurring within the first 100 ns, over the entire tem- 
perature range. This suggests that some portion of our 
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FIG. 2: (Color online) (a) ZF-/1SR time spectra collected at 
1.6 K, 50 K and 250 K for LL(C 5 H 5 N) !/ Fe2- z Se2. LF-/jSR 
time spectra collected at 1.6 K and 200 Oe. The solid lines 
are the fits to the data using the equation described in the 
text, (b) Internal field distribution, Bi„t in the vortex state 
of Liz (C5H5N)yFe2- z Se2 for the two different applied fields 
of 30 and 34 Oe. 



sample is purely magnetic in nature. To calculate the 
volume fraction of the magnetic part, we have fitted the 
following equation to the ZF data: 



A^(t) = v m A^(t) + (l-v m )A^(t) 



where 
A z J(t) = A ZF (0) 

AZ(t) = A ZF (0) 



2 1 

- exp(-A T i) + - exp(-A L <) 



(1) 



(2a) 



- + -(1 - a 2 t 2 - At) exp( — - At) 

(2b) 



Here m and nm denote the magnetic and nonmagnetic 
components, A ZF (0) is the initial asymmetry, v m is the 
magnetic volume fraction, and At, Xl are the transverse 
and longitudinal components of the muon decay for the 
magnetic fraction of the sample, respectively. A z ^(t) 
is the combination of Lorentzian and Gaussian Kubo- 
Toyabe (LGKT) relaxation function [36], [13] for the non- 
magnetic fraction of the sample with parameters a and 
A, the muon depolarization rates arise due to the con- 
centrated nuclear dipole moments and randomly orien- 
tated diluted local electronic moments, respectively. The 
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TABLE I: Parameters extracted from the fits using 
Eqs. [T] and [2] to the zero-field- /^SR data collected above and 
below T c for Li I (C 5 H 5 N) !/ Fe2- z Se2. 



T(K) 


v m (%) 


A 


1.6 


49.05±0.6 


0.18±0.01 


50 


50.8±0.6 


0.19±0.01 


250 


49.1±0.6 


0.20±0.01 



fits yield the parameters shown in Table |TJ Parameters 
obtained from the fitting show that nearly half of the 
sample is antifcrromagnetically (AFM) ordered. This is 
also consistent with other Fe-based superconductors [38| . 
The low relaxation rate A most probably stem from the 
diluted ferromagnetic Fe impurity clusters. This hypoth- 
esis is supported by the fact that it is possible to decouple 
this slow relaxation with the application of a small LF of 
200 Oe (see Fig. [2] (a)). This also proves that the dipole 
field at the muon site is static on the time scale of the 
/iSR window. The magnetic volume fraction is practi- 
cally temperature independent between 1.6 and 250 K. 
This strongly resembles the situation of phase separation 
in the alkali-metals intercalated FeSe 0, [38| . Our data 
show that there is a clear phase separation between the 
magnetic and the nonmagnetic fractions and hence the 
magnetic fraction is not affecting the properties of the 
nonmagnetic part of the sample. 

We have also performed field shift (FS) measurements 
on Lite's H 5 N) a Fe2- z Se 2 . Figure [2] (b) shows the in- 
ternal field distribution, Bi nt in the vortex state of 
Li a: (C5H5N) y Fe2-zSe2 for the two different applied fields 
of 30 and 34 Oe. FS measurements show that an appli- 
cation of additional 4 Oe magnetic field below T c does 
not shift the internal field distribution by a significant 
amount (« 0.5 G only) due to vortex pinning. This shows 
that the paramagnetic volume fraction of the sample be- 
comes fully superconducting below T c . 

To reveal the superconducting properties on 
Li a; (C5H5N)j ; Fe2- z Se2, we have performed a TF-/^SR 
study. It was very difficult to extract the supercon- 
ducting properties of Li,;(C5rl5N),,Fe2- z Se2 using a 
higher magnetic field due to higher relaxation from the 
diluted ferromagnetic impurities present in the sample. 
Therefore, we have used a small magnetic field of 30 Oe 
for this TF-/iSR study. It is noteworthy to mention 
that the lower critical field of all the iron chalcogenide 
superconductors are very low and in our case, 30 Oe 
was high enough to drive Li E (C5rI 5 N) y Fe2- z Se 2 into 
a well defined vortex state. Figure [3] (a) shows the 
temperature dependence of the internal magnetic field 
at the muon sites in Li :E (C5rl5N)j / Fe2-zSe2. The dashed 
line is drawn as a guide for the eye. We observed an 
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FIG. 3: (Color online) (a) Temperature dependence 
of the internal magnetic field at the muon sites in 
Li c (C5H5N) H Fe2-zSe2. The dashed line is simply a guide 
to the eye. (b) Temperature dependence of A -2 for 
Li;c(C5H5N) H Fe2- z Se2. The curve (black line) is a fit to the 
data using two s-wave components, each with an isotropic 
gap. 



expected diamagnetic shift of the internal magnetic field 
experienced by the muons below T c . 

To extract the superfluid density of 
Li a; (C5H5N)j ) Fe2- z Se2, we have fitted the TF muon time 
spectra using an oscillatory decaying Gaussian function, 

A TF (t) = A TF {0) exp {-a 2 t 2 / 2) cos (^B mt t + 0) 

+^(0)008(7^^ + ^), (3) 

where A TF (0) is the initial asymmetry, jg /2ir = 
135.5 MHz/T is the muon gyromagnetic ratio [30j, Bi nt 
and Bbgd are the internal and background magnetic field 
at the muon sites, is the initial phase offset of the 
muon precession signal, respectively and a is the Gaus- 
sian muon spin relaxation rate. It can be written as 
i 

a = {a 2 c + c^ m ) 2 , where a sc is the superconducting con- 
tribution to the relaxation rate due to the field variation 
across the flux line lattice and cr nm is the nuclear mag- 
netic dipolar contribution which is assumed to be con- 
stant over the temperature range of the study. 

In a superconductor with a large upper critical field 
and a hexagonal Abrikosov vortex lattice, the Gaussian 
muon-spin depolarization rate a sc is related to the pene- 
tration depth A by the expression 
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(T) $ 2 
v ; = 0.00371 



A 4 (T) ' 



(4) 



where $ = 2.068 x 10~ 15 Wb is the flux quantum [3l|]. 
Figure [3] (b) shows the temperature dependence of A~ 2 
for Li x (C 5 H 5 N) a Fe 2 _ z Se 2 . 

The temperature dependence of the penetration depth 
of Li r (C5H R N)„Fe2- j : Se 2 can be fitted using a two-gap 
s + s-wave model KM Eol , 



A" 



A- 2 (0) A- 2 (0,Ai) 



A- 2 (0,A 2 ) 



(5) 



where A (0) is the value of the penetration depth at 
T = K, Aj is the value of the i-th (i = 1 or 2) su- 
perconducting gap at T = K and u> is the weighting 
factor of the first gap . 

Each component in Eq. [5] can be expressed within the 
local London approximation 42j, |43j as 



a- 2 (r, Aj) 

A" 2 (0,A,) 



df_ 

dE 



EdE 



E 2 - Ai (Tf 



(6) 



where / = [1 + exp {E/HbT)] is the Fermi function, 
and Ai (T) = Aid (T/T c ). The temperature dependence 
of the gap is approximated by the expression 6 (T/T c ) = 

tanh|l.82 [1.018 (T c /T-l)] 51 | 

The curve shown in Fig. [3] (b) is a fit of the two- 
gap s + s-wave model to the data. The fit yields Ai — 
6.82(92) meV and A 2 = 0.93(7) meV with u = 0.30(12). 
The ratio of the larger to the smaller gap, Ai/A 2 ~ 
7.3(2) found in LL J (C5H 5 N) y Fe 2 _ z Se 2 is larger than the 
corresponding value (4.2) seen in FeSeo.85 but is consis- 
tent with the value (7.7) for the lithium amide and ammo- 
nia intercalated FeSe determined by measuring the mag- 
netic penetration depth using the [iSR, technique [U [25| ■ 
The large gap to the T c ratio is 2A 1 /k B T c = 3.96(54) 
and the small gap to the T c ratio is 2A 2 /fcsT c = 0.57(4). 
These suggest that the large gap lies above the strong- 
coupling limit whereas the small gap are in the weak- 
coupling limit. Table [TTI summarises the superconducting 
gap to the T c ratios for different iron chalcogenide super- 
conductors by means of /iSR study. 

The value of the penetration depth at T = K is found 
to be A(0) = 485(21) nm for Li^CsHsN^Fea-zSea- 
For an anisotropic polycrystalline samples, A is related 
to A a fc by A = i^Xab [45(. With this assumption, we 
found A afc (0) = 369(16) nm for Li K (C 5 H 5 N) y Fe 2 _ 2 Se 2 . 
These values are comparable with those obtained by 
Khasanov, et al. [IH for FeSeo.85 but 50 % higher than 
the value obtained from Lio.6(NH 2 )o. 2 (NH3)o.8Fe 2 Se 2 



by Lucas, et al. [21( using /iSR technique. This dis- 



TABLE II: Superconducting gap to the T c ratios for different 
iron chalcogenide superconductors by means of /iSR study. 



Compounds 


2A 1 /k B T c 


2A 2 /fc fl 


T c Refs. 


FeSeo.gs 


4.49(6) 


1.07(2) 


[25] 


FeTeo. 5 Se .5 


4.19(16) 


1.40(9) 


[44] 


K .74Fei.66Se2 


4.7(2) 




[38] 


Rbo.77Fe1.6i.Se2 


5.5(2) 




[38] 


Lio.6(NH2)o.2(NH 3 )o.f 


Fe 2 Se 2 5.27(21) 


0.69(2) 


[2JJ 


Li^CsHsN^Fea-.Se 


2 3.96(54) 


0.57(4) 


current 



crepency may be due to the fact that the interlayer 



spacing between the FeSe layers in our compound (lat- 
tice parameter, c = 23.09648 A) is 40 % higher than 
Lio.6(NH 2 )o. 2 (NH 3 )o.8Fe 2 Se 2 (c = 16.5266 A). Recent 
/iSR studies on bismuth-based high-T c cuprate supercon- 
ductors by Baker et al. [46] suggest that the penetration 
depth in this class of materials increases with increasing 
layer separation between the Cu0 2 layers with a relation 
1/d oc 1/A 2 b , where d is the interlayer distance. It is 
also found that the T c value depends only on the 2D su- 
pcrfluid density but not on the bulk supcrfluid density. 
For Li 2; (C 5 H5N) !/ Fe 2 _ z Se 2 , we obtain d/X 2 ab = 1.84(14) x 
10 4 m _1 , which is consistent with the value 2.6 x 10 4 m _1 , 
estimated for Li .6(NH 2 )o. 2 (NH 3 )o.8Fe 2 Se 2 . The nearly 
constant ratio of rf/A 2 h for the two intercalated FeSe ma- 
terials point to a noteworthy similarity to the high-T c 
cuprate superconductors. Like the Cu0 2 layer in the 
high-r c cuprates, it is the 2D FeSe layer which gener- 
ates superconductivity in this system and with an in- 
creased interlayer distance, the superfluid density within 
the layer does not change. 

In summary, /iSR measurements have been performed 
on the alkali metal-organic solvent intercalated iron se- 
lcnidc superconductor Li a; (C5tl5N)j ) Fe 2 _ z Se 2 with a T c of 
40(1) K. The superfluid density of this compound is found 
to be very low compare to other iron-chalcogenide super- 
conductors probably due to its 2D nature. Nearly 50 % 
(volume fraction) of the sample is magnetically ordered. 
A spatial phase separation has been observed between 
the superconducting and the magnetic fraction of the 
sample. The temperature dependence of the magnetic 
penetration depth is found to be compatible with a two- 
gap s+s- wave model with gap values Ai = 6.82(92) meV 
and A 2 = 0.93(7) meV. We obtain A(0) = 485(21) nm 
at T = K. Further studies are in progress to explore 
the magnetic properties and if there is any microscopic 
coexistence of superconductivity and magnetism in the 
magnetic volume of this material. 

The /zSR experiments were performed at the Swiss 
Muon Source, Paul Scherrer Institut, Villigen, Switzer- 
land. A.K-M gratefully acknowledges the financial sup- 
port from Institute of Physics University of Zurich. 
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